Potential and current distribution in strongly anisotropic Bi 2 Sr 2 CaCu 2 08 single 

crystals at current breakdown 
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Experiments on potential differences in the low-temperature vortex solid phase of monocrystalline 
platelets of superconducting Bi2Sr2CaCu2 08 subjected to currents driven either through an ab 
surface or from one such surface to another show evidence of a resistive/nonresistive front moving 
progressively out from the current contacts as the current increases. The depth of the resistive region 
has been measured by an in-depth voltage probe contact. The position of the front associated with an 
injection point appears to depend only on the current magnitude and not on its withdrawal point. It 
is argued that enhanced nonresistive superconducting anisotropy limits current penetration to depths 
less than the London length and results in a flat rectangular resistive region with simultaneous ab 
and c current breakdown which moves progressively out from the injection point with increasing 
current. Measurements in ab or c configurations are seen to give the same information, involving 
both afo-plane and c-axis conduction properties. 

PACS numbers: 74.72.Hs, 74.25. Fy, 74.25.Sv 



I. INTRODUCTION 

Understanding the current and the potential distri- 
butions at the onset of dissipation in a superconduc- 
tor is important; at a fundamental level for interpreting 
the results of transport measurements in terms of the 
force- velocity relation for the vortices and the Josephson 
coupling between layers as well as on a practical level 
for understanding how to maximize the current carry- 
ing capacity of a wire. The family of high-T c cuprate 
superconductors naturally raises the question of the in- 
fluence of high anisotropy on current carrying proper- 
ties. How does one interpolate between the extreme 
anisotropic limit, where current injected into a super- 
conducting plane has no transfer to other planes, and 
the isotropic case where the penetration of the current in 
the superconducting state is determined by the London 
screening lengthi To give a concrete illustration of the 
problem, a naive interpretation of the threshold current 
for dissipation in a strongly anisotropic Bi2Sr2CaCu20s 
(BSCCO) single crystal based on penetration limited to 
the London screening length or, even more naively, sup- 
posed to be uniform over the typical thickness of a few 
micrometers gives threshold current densities one to two 
orders of magnitude lower than a standard interpretation 
of a magnetic hysteresis loop. 

Anisotropy is a central feature of the family of the 
high-T c cuprate superconductors which are usually mod- 
eled by discrete superconducting sheets parallel to the 
ab crystallographic plane weakly coupled together in the 
c direction by extended Josephson junctions? 2 " The weak 
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coupling between planes gives rise to very high conduc- 
tion anisotropy in both the normal and superconduct- 
ing states. It is also responsible for the richness of the 
phase diagram of the vortices created when one applies 
a magnetic field. A convenient experimental probe for 
investigating these phases, particularly as regards their 
pinning to the host lattice disorder, is to look at the volt- 
age response to transport current (VI characteristics) for 
currents which take one into the resistive regime where 
the vortices are dislodged. But to interpret these mea- 
surements correctly, it is important to know about the 
current distribution in the resistanceless regime, how it 
is modified as dissipation sets in and of course what the 
potential distribution then looks like. 

It was remarked some time ago* 3 - in connection with 
resistivity experiments, that if the response of the super- 
conductor is Ohmic, E = p • J (the resistivity tensor p 
is constant and independent of current and position), the 
problem may be treated as for the normal state where it 
is implicit that the phase slip giving rise to the voltage 
also relaxes the magnetic dephasing to allow the current 
to penetrate beyond the London screening depth. In this 
regime, at least, the high anisotropy of the (BSCCO) 
samples studied limited the current penetration into the 
depth of the sample to a small fraction of its thickness 
and that it would therefore be quite erroneous to deduce 
a critical current density on the basis of uniform penetra- 
tion across the section of a short sample. Experiments 
to date on the current profile have been able to show a 
concentration of current towards the sample edges, but 
have been insensitive to the depth dependence^ which is 
the aspect that interests us chiefly here. 

The present paper reports an experimental investiga- 
tion and a few basic reflections on this problem in the 
non-Ohmic, low-temperature, high magnetic-field vor- 
tex solid phase, also in monocrystalline BSCCO in a 
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FIG. 1: a6 and c contact configurations. 



c-directed magnetic field. The experiment consists of 
measuring the voltage response to short pulses of cur- 
rent up to and beyond dissipative breakdown. Potentials 
are measured at contacts placed in the usual way on the 
ab faces supplemented by in-depth potential contacts to 
have direct access to depth-dependent features. One of 
the important points that emerges is that the physically 
unrelated c-axis and a6-plane dissipative breakdowns in- 
teract to create a resistive/nonresistive front which has 
the effect of altering the current distribution and giv- 
ing rise to a critical current which involves both ab and 
c characteristics. Not only does the voltage response 
along an afo-plane surface into which the current is in- 
jected and withdrawn (the "ab configuration" of Fig. 
show features of the Josephson junctions between planes, 
but also the response to current injection and withdrawal 
on opposing a6-plane faces (the "c configuration" of Fig. 
^) shows features of breakdown in the ab plane. One 
manifestation of this is that the threshold current in the 
ab configuration on BSCCO single crystals 5 shows the 
same temperature and preparation dependence as for the 
c configuration^ This brings up the question as to who 
measures what and does one measure what one thinks; 
but beyond that what do the results mean? 

II. EXPERIMENTAL TECHNIQUE 

Our samples are thin monocrystalline platelets with 
the c axis perpendicular to the face. As all the sam- 
ples studied gave qualitatively the same results for simi- 
lar contact configurations, we limit ourselves here to the 
results on one sample of each different contact config- 
uration. The sizes of samples A, B, and C are about 
2 x 0.35 x 0.005 mm 3 , 1.1 x 0.3 x 0.002 mm 3 , and 
0.7 x 0.5 x 0.005 mm 3 . They were all fabricated by a melt 
cooling techniquei^ The critical temperatures were be- 
tween 88 and 90 K with a transition width of about 2 K at 
zero field. The anisotropy coefficient 7„ w 500 was esti- 
mated from normal-state resistivities with p a b ~ 100 fj,Q 
cm at 90 K. In the case of samples A and B, electrical 
contact was made by bonding 25-fim gold wires with sil- 
ver epoxy fired at 900 K. For sample C the contacts were 
made by depositing silver on a lithographically defined 
area using a lift-off technique and heat treating at 700 K 
for 1000 s. In all cases the contact resistance was less 
than 3 £1. The sample was mounted flat against a 7-mm- 
diameter 0.5-mm-thick sapphire disk with silicone grease 
to ensure thermal homogeneity and mechanical freedom 



and placed in the bore of a superconducting magnet with 
the field along the c direction. The sample and disk 
were surrounded by exchange gas and the temperature 
was electronically regulated. The longitudinal voltage- 
current (VI) characteristics were measured by a sym- 
metrized differential four-point technique using 25-fj.s tri- 
angular pulses (12.5 fj,s from zero to maximum current) of 
maximum amplitude in the range 10 mA < I m < 350 mA, 
usually restricted to exceed the threshold by about 30%, 
with repetition time 0.1 s. As reported earlier^ the 
observed threshold is independent of pulse duration for 
times up to about 250 /is, argued there to be the time 
required for heat produced in the current contact to dif- 
fuse towards the voltage contact, whereas bulk heating in 
the superconductor was estimated not to be important. 
Further technical details may be found in Ref. Oj. 



III. EXPERIMENTAL RESULTS 

We present first the results of experiments on the 
VI response with surface contacts in the standard ab 
and c configurations on the same sample using con- 
tacts common to the two configurations. The results 
lead rather naturally to a description in terms of resis- 
tive/nonresistive fronts moving outwards from the cur- 
rent contacts. Experiments are then presented for other 
surface contact configurations intended to check this hy- 
pothesis. These are followed by experiments with a dif- 
ferent type of contact configuration designed to probe the 
depth of the sample as breakdown progresses. 

The first experiment was done on samples contacted 
to be able to measure in both ab and c configurations 
on the same sample at the same time as illustrated in 
Fig. [3 The "ab configuration" measurements refer to 
the potential drop between contacts (A2, A4) for current 
injected at Al and withdrawn from A5. The "c configu- 
ration" refers to the potential measured across (A2, Al) 
for current passed through (A1,A&). Figure shows the 
temperature dependence of the threshold current in the 
VI characteristic at 1.5 T in sample A for both configu- 
rations. The threshold current is defined by the break in 
slope criterion of Ref. which gives essentially the same 
values as the criterion of crossover from Kim- Anderson 
behavior at low current to power law close to linear be- 
havior at higher current, 9 illustrated in detail for the first 
pair of curves of Fig.|5|in the semilog plot in the same fig- 
ure. The measured thresholds, as represented in Fig. |3 
are seen to be independent of the configuration, whether 
the sample is field cooled (FC) or zero-field cooled (ZFC) 
prepared, despite the fact that the two preparations give 
different thresholds for temperatures lower than the peak 
in the ZFC results^ For the FC preparation, the field is 
applied above T c and the sample is cooled at constant 
field to the lowest measuring temperature, subsequent 
measurements being made on increasing the tempera- 
ture at the same field. For ZFC preparation, the sys- 
tem is cooled in zero field from above T c to the lowest 
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FIG. 3: Temperature dependence of the threshold current 
with ab and c contact configurations on the same sample. 
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FIG. 2: VI curves measured on sample A. The letters on the 
VI curves refer to the contacts used: the first pair in brackets 
denote the current contacts and the last two the potential 
contacts. The inset in (a) shows the contact configuration. 
Part (a) shows the response to up and down current sweeps of 
isoceles triangular shape of total duration 25 /is. Parts (b) and 
(c) show in more detail, on both linear and semilogarithmic 
plots, upsweep response for the first two curves of part (a). 
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FIG. 4: Magnetic-field dependence of the threshold current 
with both contact configurations on the same sample. The 
solid line shows H^ 1 ^ 2 . 



temperature before applying the field and then making 
measurements at sequentially higher temperatures at the 
same field. Figure^shows the magnetic-field dependence 
in both contact configurations for ZFC preparation at 
5 K on increasing the field (FC preparation followed by 
field variation gives the same result, except for the ini- 
tial FC prepared pointjS). At fields above about 0.3 T, 
the threshold currents for the two different configurations 
again show the same behavior, this time in field, varying 
approximately as if -1 / 2 . Not only are the temperature 
and field dependences of the threshold currents in the two 
configurations similar in form, but the values themselves 
coincide. Measurements along several other lines in the 
(H, T) plane confirm this indistinguishability of ab and c 
configuration thresholds as generic behavior for the low- 
temperature high-field domain. This fact led us to look 
in more detail with the following experiment. 

Again referring to the inset of Fig. |2 we used three 
potential contacts (A2, A3, A4) on the top layer and two 
(A7, A8) on the bottom. The contact A2 is half as far 
away from the left current contact ^41 as is the contact 
AA from the right-hand current contact A5, while AS is 
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FIG. 5: VI response of sample B. The labeling convention is 
as before. 
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FIG. 6: VI response for sample C. The dimensions k and h 
are 125 /jm and 75 /im, respectively. The potential contact 
dimension along the sample length is 50 /im. 



midway between the two current contacts. If we apply 
the current between contacts Al and Ab (ab direction) 
we measure a much lower threshold current between A2 
and A4 than between ^43 and AA. On the other hand, 
we find virtually identical VI curves on the ab potential 
contacts (A2,A4) whether the current is applied in the 
c direction through contacts (A1,A6) or along the ab 
direction (Al,A5); similarly the VI across (A3,A4) is 
nearly identical for afe-directed current through (Al,A5) 
or for c-directed current through (AB, A9). 

These results suggest the idea that the sample first 
becomes resistive around the injection points and that 
this resistive region progressively invades the sample as 
the current is increased. This led us to experiment with 
the contact configuration of sample B illustrated in Fig. 
where we drove the current through the c direction 
and also measured the voltage drop across that direction. 
On sending the current through (B1,B4) the VI char- 
acteristic measured across (B2, B5) shows a threshold 
current which is about half that measured across con- 
tacts (B3,BQ). On the other hand, if we put the cur- 
rent through the middle contacts (B2, B5) we find the 
same threshold at either of the two end potential contacts 
(Bl, B4) or (B3, B6) with a value intermediate between 
the previous two. 

Although these experiments lend considerable support 
to the idea of a resistive/nonresistive front progressing 
along the surface, they say nothing about its profile 
with depth since at no point were we able, in this low- 
temperature high-field phase and despite micron thin 
samples, to drive the lower face contacts resistive with 
current sent along the top face. To have more direct in- 
formation about the shape of the front with depth, we 
prepared a third sample C as in Fig. with a lithograph- 
ically defined terrace argon ion etched along the length 
of one ab face to a depth of 220 nm. The upper part of 
the face was also argon ion etched to avoid any difference 
in surface pinning between it and the terrace, the width 
of which was restricted to a small fraction (20%) of the 



total sample width to minimize the trivial effect of cur- 
rent spread as the distribution attains the depth of the 
terrace. The contacts were accurately aligned between 
the two levels by a scanning electron-beam masking mi- 
croscope. On a control sample with the same contact 
configuration but with no terrace, the VI response on 
contacts (C2, C3) gave the same values of the thresh- 
old current and dynamic resistance as contacts (C5, C6), 
themselves very similar to those for the top contacts of 
the terraced sample. 

Sending an afe-directed current through contacts 
(Cl, C4) on the upper part of the face, we measured the 
potential drop across (C2, C3) on the upper part and 
across (C5, C6) on the step. The threshold current for 
the terrace contacts is considerably higher than for the 
corresponding top contacts: about 25 mA and 10 mA, 
respectively, much greater than the decrease of 20% in 
the current density to be expected from simple current 
spread into the wider part of the sample. We interpret 
the large difference to arise from the current distribution 
with depth. 



IV. INTERPRETATION OF EXPERIMENTAL 
RESULTS 

All these observations are consistent with the idea of 
a resistive front moving progressively outward from the 
current contacts as illustrated for sample A in Fig. [7| 
Also, and more surprisingly, it seems that the position 
of the front depends only on the magnitude of the cur- 
rent injected and seems to be independent of its ultimate 
destination if one is to account for the nearly identical 
VI response on, for example, the potential contact pair 
(A2, AA) independently of whether the current is with- 
drawn on the same face (A1,A5) or the opposite face 
(A1,A6); similarly for the potential pair (A3,AA) for 
current through (A1,A5) or (A5,A9). The difference 
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FIG. 7: Schematic illustration of a lengthwise section of sam- 
ple A showing the propagation of resistive fronts with current 
directed in the (a) ab direction (Al,A5) and (b) c direction 
(Al,AQ) or (A5,A9). The c axis is vertical. 
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FIG. 8: Schematic illustration of the propagation of resistive 
fronts with a c-directed current in a lengthwise section of sam- 
ple B. The current was driven through (a) (Bl, B4) contacts 
near one end of the sample (b) (B2, B5) at the middle of the 
sample. 
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FIG. 9: Schematic illustration of a section of sample C with 
an etched terrace to investigate depth dependence of the re- 
sistive/nonresistive front which propagates from current con- 
tact CI to attain the top surface contact C2 and, at higher 
current, the terrace contact C5. If the front were rectangu- 
lar as illustrated and had a current-independent aspect ratio 
£/d, the dimensions would be linear in the total current and 
£/d= (VdiX^/'-O provided that 1% > 1%. 



to progress in the experiment with sample B when the 
current is run through the c direction. Here too the 
lower threshold current was obtained for potential pair 
(B2, Bb) nearest the current contacts. The resistive front 
from current contacts (SI, BA) attains (B2,Bb) before 
(B3, B6). When the current was run through the middle 
of the sample by (B2, Bb) the injected current is divided 
between two directions with a corresponding reduction 
in current density so that the voltage drop appears at 
higher current for the same distance and gives the same 
response at both ends. 

Figure [5] represents the front for the experiment with 
the terraced sample. Here the form of the front has been 
taken to be rectangular for the reasons outlined in the 
following section. As the front is defined by the line at 
which the current density reaches the threshold for re- 
sistive breakdown (ab plane and/or c direction) the inte- 
gral of the threshold current density crossing the front is 
just the total current. Furthermore if the front does not 
change geometrical form with current amplitude its di- 
mensions increase linearly with the current. Then, know- 
ing the distance from the nearest current contact and the 
depth of the terrace, the threshold current values for the 
upper and step contacts give information on the depth 
and potentially some information on the form of the front 
as it attains the upper and terrace level contacts. We can 
expect the aspect ratio of the front to be related to the 
anisotropy. 



between the values of the threshold current measured 
on (A2, AA) and (A3, AA) contact pairs is accounted for 
by the distances of the potential contact from the near- 
est current contact (A2 from Al or AA from Ab). In 
these experiments, the fronts did not attain contact A3 
which remained at the potential of the resistanceless re- 
gion. The triangular form used for the resistive region 
in the figure is purely schematic to show the progression 
of the front. Further consideration of possible forms is 
given below where it is argued that the form is in fact 
most probably rectangular. 

Figure El represents how the resistive front is imagined 



V. ELEMENTARY UNDERSTANDING 

We shall assume for the present discussion that edge 
effects are not a dominant factor for our observations. 
Although we expect an integrable singularity at the edge 
of the sample, a large part of the middle of the sam- 
ple has an essentially flat distribution with distance from 
the center line^S and we deal in the first instance uniquely 
with the distribution with depth z and along the length x 
of the sample. We know that in the normal resistive state 
the current distribution is governed by the anisotropy of 
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the resistivity. By rescaling t he coor dinates according to 
the anisotropy factor j n — y/ p c / p a b the problem can be 
solved as a Laplace equation in the potential. The neces- 
sary coordinate rescaling is determined by the anisotropy 
factor 7: z — > 7Z. In the normal state, it is clearly given 
by 7 2 = p c / Pab ~ (500) 2 and would give rise to an ab 
current penetration depth of about p a b ~ £ a b /in ~ 1 /im 
for a distance £ a b — 1 nrm between ab current contacts^ 

It is less recognized, perhaps, that the low current dis- 
tribution in the nonresistive superconducting regime can 
be governed by the same sort of equation for the super- 
conductor phase in similarly rescaled coordinates if the 
penetration depth is smaller than the London screening 
depth. In the nonresistive superconducting phase, the 
Laplacian in the potential is replaced by a Laplacian 
in the superconductor phase which has the same form 
of relationship to the superconducting current as does 
the potential to the current in the normal state. The 
anisotropy ratio in the superconducting state is given by 
7! = (n s eh/m a b)/sJj, the ratio of the coefficients of the 
ab and c phase gradients in the relation between the cur- 
rent operator and the phase. The sJj term is the low 
current expansion of the Josephson relation between suc- 
cessive superconducting planes separated by s, while n s 
is the three-dimensional (3D) superconducting electron 
density and 1/m the inverse mass tensor. If one were to 
neglect London screening, 7 S would govern the current 
penetration in the same way as does 7„ in the normal 
state. 

The superconducting anisotropy is usually expressed as 
a mass ratio or London screening length ratio for c and 
ab currents 7 2 = m c /m a b = A 2 /A 2 6 and is customarily 
expected to have about the same value as the extrapo- 
lated normal-state resistivity ratio. In BSCCO, however, 
experiments 1 ^ on mesa structures have shown a Joseph- 
son critical current reduced by a factor of about 30 with 
respect to the usual Ambegaokar-Baratoff 11 relation to 
the high current resistance; the low current resistive slope 
is proportionately altered to V gap /sJj ~ 30p c and the 
coefficient in the Josephson relation J z — Jj sin(</J n — 
fn+i) ~^ sJjdip/dz at low current is reduced by the 
same factor with a corresponding enhancement of the 
anisotropy factor [tp n is the phase ip(r) of the supercon- 
ductor wave function at the nth superconducting plane] . 
This has been attributed to the d-wave nature of the or- 
der parameter and already enhances the anisotropy factor 
by V30 over the normal state. But if one also considers 
the effect of phase fluctuations across the Josephson junc- 
tion caused by misalignment of vortices in magnetic field, 
another multiplicative factor (cos ^n.n+i} -1 ' 2 should be 
incorporated into 7 S . The factor (cos(^ n! „+i) has been 
variously estimated from Josephson plasma resonance 12 
experiments to be between about 10 -3 and 10 -1 , ex- 
trapolating to about 4 x 10~ 3 (Refs. and 0) at low 
temperature and poH = 1.5 T, which suggests a very 
strongly enhanced anisotropy parameter in the super- 
conducting phase as compared with the normal state: 
3 x 10 4 < 7 S < 1 x 10 6 . Current penetration in the super- 



conducting phase is governed by a combination of Lon- 
don screening and anisotropy, both of which act to limit 
the penetration. But, even before including the effects of 
London screening, the anisotropy limits the penetration 
depth p a h ~ tab/ Is in the middle of a 1-mm-long BSCCO 
sample to 2 < p a b < 20 nm <C X a b — 200 nm. Thus in 
short samples of BSCCO like ours the current penetra- 
tion in the superconducting phase is always dominated 
by anisotropy rather than by London screening. 

This suggests that the effective thickness for current 
flow might be given approximately by p a b ~ P-ab/ls < 
A a 6 <C t (the real thickness), and furthermore that the 
current density is nonuniform along the sample, being 
higher closer to the contacts, causing breakdown to re- 
sistive behavior to be progressive and not to penetrate 
immediately throughout the bulk. Under these condi- 
tions, ab breakdown alone is not possible: if only ab 
breakdown were to occur, the front would have to be 
parallel to the c direction because no part of the inter- 
face with the nonresistive region can be parallel to the 
electric field produced along the ab breakdown direction. 
But if breakdown does not occur uniformly throughout 
the thickness of the sample, the front must somewhere 
change orientation. The electric field in the resistive por- 
tion, which must be oriented perpendicular to the inter- 
face, therefore also changes orientation and that requires 
c breakdown. Hence a resistive/nonresistive front within 
the sample which does not go straight through it along 
a principal direction necessarily implies simultaneous ab 
and c breakdown. The regions of ab and c breakdown 
need not be spatially coincident however, but they must 
be at least contiguous. The same can be argued for the 
case of pure c breakdown. 

Another consequence of the nonuniformity of the cur- 
rent distribution in depth is the presence of shear forces 
on the vortices and the possibility of sliding between 
planes of (semi)ordered vortex segmentsA 1 ^ This has 
important consequences on the locality of the E(J) 
relation. If the force required to shear two adjoin- 
ing planes derives from an effective Josephson coupling 
Jj(cos y>n,n+i)) the condition for nonelastic shear (slip) is 
dJab/dz > Jj(cos ip n<n +x)a/ s 2 , where a and s are the dis- 
tances between vortices in the plane and the distance be- 
tween neighboring planes, respectively. Using the expres- 
sion for the anisotropy penetration depth p a b ~ £ a b/"fs 
and identifying the effective Josephson coupling with the 
c-axis critical current density J* h , the condition for shear 
to occur at the onset of dissipation is J^/J^ > jalab/s 2 . 
At high values of anisotropy, j s > 10 2 , however, the 
shear strength is dominated by magnetic coupling be- 
tween vortex segments i&lL Either by relating the crit- 
ical force per vortex segment for shear slip f c to the 
tilt modulus C44 by f c » C^as and using the evalua- 
tions of the modulus for magnetic coupling^ or using 
directly the calculation of the shear strength,- 8 it ap- 
pears that shear slip between planes should occur for 
dJ a b/dz > (aj)o/32n 3 \* b )(cos<p n , n +i)a/s. It will be seen 
that this condition is met with the values obtained from 
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the data analyzed according to the reasoning outlined be- 
low and it is a necessary condition to have a local E(J) 
relation. 

The above considerations lead to a scenario where the 
current injected and withdrawn from the same ab surface 
is confined by anisotropy to a penetration depth ~ £ ab /"f s 
in the nonresistive regime with consequent enhancement 
of the current density. When the current density attains 
the critical value, for either J c or J ab , the current distri- 
bution is modified until both reach their breakdown val- 
ues and the sample is resistive on one side and nonresis- 
tive on the other side of a front which advances out from 
each current contact point as the current is increased. 
As the fronts attain successive voltage contacts, these 
display a potential drop with respect to the portion of 
the sample which has remained in the nonresistive state. 
The shape of the resistive region can be expected to be 
a function of the E-J response around the breakdown 
values (vortex depinning for ab current and Josephson- 
junction behavior for c current) described for the most 
part by the anisotropy factor 7 of the resistive state and 
r = J*5>/^c h - If tne shape were determined uniquely by 
the E-J response, its position would depend only on the 
contact from which it originated and the magnitude of 
the current injected. 

Restrictions are imposed on the shape of the front by 
the irrotational nature of the electrochemical field and 
by current conservation. The inductive field from the 
current is negligible compared to resistive fields for our 
pulses and sample sizes and we can suppose the field E 
to be irrotational: VxE = — dti/dt = 0. We can then, 
as in London's treatment of breakdown in a type-I cur- 
rent carrying wireji demand that § E-d£ = be satisfied 
through the resistive/nonresistive interface fFig. HUI) with 
the consequence that the interface must be perpendicu- 
lar to the field. As argued above, our quasipoint contact 
situation enforces that both ab and c breakdown appear 
simultaneously and the VI response in either configura- 
tion will contain both c and ab features. Either ab and c 
breakdown occur together at the interface or they occur 
alone in separated but necessarily contiguous or overlap- 
ping regions. The first case is considered in part (a) of 
Fig. ^5] where a segment of the interface is locally ori- 
ented at an angle 8 with the ab plane. A model with a 
local (locality implicitly assumes shearing of the vortex 
segments from plane to plane) E-J relation for J > J th 
of the form 



E x = 0, 



\J,\ < J. 
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FIG. 10: The upper part (a) of the figure illustrates the condi- 
tion that the electric field be perpendicular to any segment of 
the resistive/nonresistive interface supposing spatially coinci- 
dent ab and c current breakdown; "res" denotes the resistive 
side, "nonres" the nonresistive side. E a t, and E c are given 
by the E(J) relationship, an example of which is shown in 
part (b). Positive differential resistivities for both compo- 
nents imply both dJ a t/dx a i, > and dJ c /dx c > 0, in contra- 
diction with current conservation which requires that these 
two derivatives have opposite signs (for a 2D distribution). 
If however one of the two current-density components is sub- 
critical, O = 0, it/2, the corresponding derivative is free. The 
latter situation is illustrated in the lower part of the figure 
where in addition the front changes orientation from O = 
to 6 = 7r/2 within the sample, a situation which is still com- 
patible with the irrotationality condition <f E ■ d£ — on the 
electric field provided that the resistive region between pure 
ab and pure c breakdown also terminates at the corner of the 
rectangle. 



E x = ±Ef b + p' ab (J x T 4b), J* % ±4b (1) 
E z = 0, |J z |<J c th 

E z = ±Ef + p' c (J z T Jf), J z ^ ±Jf (2) 



on the resistive side is shown in part (b) of the figure. p' ab 
and p' c are the dynamic resistivities in the ab and c di- 
rection in the superconductive state beyond breakdown. 
The condition that the field be perpendicular to the inter- 
face would demand that tan 6 = Az/Ax = ±E t J b 1 /E t J !1 , 
where the E^ and E^ h components of the electrochemical 



field both point towards the interface, in the directions 
of the ab and c components of the current flow. However 
because both components of the current density must di- 
minish on approaching the interface, dJ x /dx and dJ y /dy 
both have the same sign whereas current conservation 
V • J = requires them to have opposite signs. More 
explicitly, if the interface is imagined to have positive 
intercepts with the axes of an (x, z) coordinate system 
whose origin is at P in the resistive region as in Fig. 1101 
both J x > J*]- and J z > J* h (all with positive values, 
current flowing from left to right and top to bottom). 
But the fact that the current diminishes on approach- 
ing the interface along each of the principal directions 
to reach J*£ and J' h which defines the interface imposes 
that both dJ x /dx < and d,J z /dz < at the same point 
P, in violation of current conservation which demands 
that dJ x /dx = —dJ z /dz. 

We conclude that a region of spatially coincident ab 
and c breakdown may not have an interface with the 
nonresistive region. 

Resistive/nonresistive interfaces perpendicular to pure 
ab or pure c breakdown, however, do not violate current 
conservation. For sole c breakdown, J z > J* h , J x < J*£ 
at the interface and the positive dJ x /dx which must re- 
sult from the negative dJ z /dz to conserve current has 
no consequence for the electric field because the x com- 
ponent is subcritical and produces no electric field. A 
similar argument can be made for pure ab breakdown, 
establishing that nonresistive/resistive interfaces are pos- 
sible for single-component breakdown whereas they are 
not for spatially coincident two-component breakdown. 
Furthermore it is possible to change the orientation of the 
interface between the two perpendicular principal direc- 
tions of the superconductor on the condition that the two 
separated mutually perpendicular ab and c fronts meet at 
a common point and that the interface between the two 
types of breakdown converges to the same point. The 
latter interface could be either a simple common bound- 
ary or more generally a lens shaped region of coincident 
ab and c breakdown. Such a situation is illustrated in 
part (c) of Fig. ^| We conclude that the fronts must be 
formed of segments of interfaces perpendicular to pure ab 
or pure c current breakdown. Numerical solutions of the 
modeliS indicate a simple rectangle with a generic form 
of the type illustrated in Fig. ^] 

The length £ and depth d of the resistive region are 
related to the current by 

I = w{djf b +£jf), (3) 

where w is the width of the sample. The threshold cur- 
rent measured for the potential between the nonresistive 
region and a contact on the surface at a distance £\ from 
the current injection contact is given by the condition 
that the front arrives at the contact, I = l\: 

r£ = wii{J%d/t + j?) 

and for a voltage contact at depth d\ by d = d\ : 
l£ = wd 1 (l/d)(J%d/£+0, 
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FIG. 11: Generic form of resistive regions. The dotted cur- 
rent withdrawal arrows are ab and c configuration alterna- 
tives. Horizontal lines indicate ab breakdown, vertical lines c 
breakdown. 



provided that £(I) is sufficient that the front attains 
the contact in the x direction (1^ > I^t). In general, 
£ = £(I) and d = d(I), but if there are no other relevant 
lengths in the problem (distances between ab configura- 
tion current contacts £ ab S> £ and c configuration current 
contacts £ c = t ^> d) one can expect the ratio £/d to 
be independent of the current and only a function of the 
anisotropy factors 7 and T = J^/J* 11 . If we interpret 
the top and terrace thresholds of the experiment on sam- 
ple C as measurements of 1^ and and we make the 
hypothesis that the aspect ratio £/d is independent of 
current, then the measurements tell us that 

*/d=(Vd!)(j£/j£)« 2500 

and 

J^d/t + Jf w 10 A cm" 2 (4) 

or, multiplying by £/d, 

Jab + Jft/ d ~ 2 -5 x 10 4 A cm" 2 , (5) 

where the last two quantities arc the same combination 
of ab and c transport properties which we cannot sepa- 
rate without a model calculation for the aspect ratio of 
the resistive front. This will be the subject of a future 
paper— on numerical solutions of the present model. 

The upper limits on J* h and J*£ given by relations Q 
and (JSJ) are consistent with the value J* h w 2 A cm -2 
obtained from combining Jj measured in the mesa ex- 
periment with (cos (pn^n+i) «4x 10~ 3 extrapolated from 
the Josephson plasma resonance experiments 1 ^ and the 
value J*£ re 2 x 10 4 A cm estimated from the magnetic 
hysteresis— The condition for slip between planes is well 
satisfied for these values indicating that the description 
is self-consistent. 
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VI. SUMMARY 

The experimental results on VI characteristics for 
ab surface contacted samples of the highly anisotropic 
BSCCO superconductor in the low-temperature vortex 
lattice phase above about 2000 Oe show that the VI 
responses for ab or c current configurations are virtually 
identical and that the dissipative region invades the sam- 
ples in the form of a resistive/nonresistive front moving 
out from the current contacts with increasing current. 
The VI characteristic depends only on current magni- 
tude, voltage contact position and nearest current injec- 
tion point. It is argued that current penetration in the 
nonresistive regime is limited by the high anisotropy to 
depths less than the London screening length, and that 
the resistive region must involve simultaneous ab and c 
breakdown and shear of successive planes of vortex seg- 
ments. The shape of the resistive region is argued to 
be composed of rectangular segments and seems to be a 
simple rectangle with interfaces to the nonresistive por- 
tion consisting of pure ab or pure c breakdown. A differ- 
ent contact arrangement using an ion etched terrace to 



sample the potential in the depth of the sample brings 
extra information that allows access to the aspect ratio 
of the rectangle. But because the ab and c configuration 
experiments measure the same combination of ab and c 
properties, it is necessary to propose and solve a specific 
model to separate them. A detailed discussion of numer- 
ical solutions of the model has been reserved for future 
publication. 
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Fuchs and collaborators (Ref. |4|) have done a very inter- 
esting experiment to determine the distribution of the c- 
directed component of the current-induced magnetic field 
across and just outside a current carrying sample. In the 
low-temperature phase which interests us here, they de- 
duce that there is a stronger current density on the edge 
of the sample. Their observations are compatible with the 
expected finite width correction j x = [1 — (2y /w) 2 ] -1 ^ 2 , 
where w represents the full width and y is measured from 
the midline. This distribution is constructed to ensure the 
zero c-directed induced field in the sample as required by 
the superconducting state with fixed vortices. 



